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The A-scaling method has been applied to ultra-relativistic p+p, C+C and Pb+Pb collision data 
simulated using a high energy Monte Carlo package, LUCIAE 3.0. The A-scaling is found to be valid 
for some physical variables, such as charged particle multiplicity, strange particle multiplicity and 
number of binary nucleon-nucleon collisions from these simulated nucleus-nucleus collisions over an 
extended energy ranging from Ei a b = 20 to 200 A GeV. In addition we derived information entropy 
from the multiplicity distribution as a function of beam energy for these collisions. 

PACS numbers: 12.38.Mh, 24.60.Ky, 24.85. +p, 25.75.Nq 



Ultra-relativistic heavy ion collisions provide a unique 
means to search for a new state of matter, the Quark- 
Gluon Plasma (QGP), where quarks and gluons over an 
extended volume are de-confined [1-3] and a phase tran- 
sition between hadrons and the QGP will occur [4-6]. 
There could be a discontinuity in global features, in par- 
ticular, the entropy from particle multiplicities, in nu- 
clear collisions associated with the onset of the phase 
transition. In this Letter we report the first application 
of the A-scaling method and information entropy calcu- 
lation in ultra-relativistic heavy ion collisions. 

In intermediate energy heavy ion collisions, A-scaling 
was proposed by Botet and Ploszajczak [7]. They applied 
the A-scaling law to the INDRA data in intermediate en- 
ergy heavy ion collisions (Xe+Sn, 25-100 MeV/nucleon) 
by using Z max (the maximum of charge in reactions) 
as order parameter [8] and found that the distribu- 
tions of Z max obey the A = 1/2 scaling law below 32 
MeV/nucleon while they obey the A — 1 scaling law 
above 32 MeV/nucleon. This indicates that a transi- 
tion from an order phase to the maximum fluctuation 
phase (disorder phase) occurs around 32 MeV/nucleon. 
Recently, Ma et al. analysed the multi-fragmentation 
data for mass around ~ 36 light nuclei systems. They 
found a similar change of A-scaling occurs when the ex- 
citation energy of the system around 5.6 MeV/nucleon. 
Combination of analysis with other probes of phase tran- 
sition (eg., Zipf law [9] and maximum fluctuations [10]) 
reflects a transition of matter phases with the change of 
A-scaling [11]. Since A-scaling is a useful tool to identify 
the phase transition, we try to make the similar analysis 
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in ultra-relativistic heavy-ion collisions. 

Botet and Ploszajczak proposed A-scaling to identify 
the transition in intermediate energy heavy ion collisions 
[12,13]. The A-scaling law is observed when two or more 
probability distributions P[m] of the stochastic observ- 
able to collapse onto a single scaling curve $(z) if a new 
scaling observable is defined by: 



This curve is: 
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where A is a scaling parameter, to* is the most probable 
value of to, and (to) is the mean of to. When A = 1, 
this kind of scaling law is called the first scaling law that 
is caused by self-similarity of system. This self-similarity 
means that, if these distributions with different (to) by a 
new kind of variable z, they entirely collapse on the same 
curve. In fact, the famous KNO scaling [14] is the special 
case that the A=l scaling law holds with a stochastic 
observable of multiplicity of particles. The INDRA data 
was explored using the A — 1/2 and 1 scaling laws and 
a phase transition was observed [15]. If we assume that 
P[m] is a Gaussian distribution, we have: 
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where /i = (m) = to* , a is the width of the Gaussian 
distribution, they both depend on incident energy. If 
this Gaussian distribution P[m] obeys A-scaling law , 
we should have: 
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On the other hand, the information entropy is observ- 
able such that it was proposed to describe the fluctuation 
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and disorder of a system by Shannon [16]. 
by: 
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It is denned 



(5) 



where J P(m)dm = 1. This indicates that the quantity 
of fluctuation of system and depend on the distribution 
P[m]. If P[m] is the Gaussian distribution, we have : 



H = Ina + (1 + ln2w) /2 ~ Ina + 1.419 



(6) 



In this work, we investigate the A-scaling law and the 
information entropy for p+p, C+C and Pb+Pb in high- 
energy collisions (20-200 AGeV) with help of LUCIAE 
3.0 of Sa and Tai [17]. The head-on collisions are sim- 
ulated in this work. The LUCIAE is of a Monte Carlo 
model and is an extension of the FPJTIOF [18]. Here 
the nucleus - nucleus collision is described as the sum 
of nucleon-nucleon collisions. LUCIAE was improved 
in the following three aspects: (1) Re-scattering of fi- 
nal hadrons, spectator and participant nucleons are con- 
sidered in the LUCIAE, because the role of re-scattering 
[19] can not be neglected in high-energy domain. (2) The 
LUCIAE implements the Firecracker Model that includes 
collective multi-gluon emission from the color fields of in- 
teracting strings in the early stage of relativistic ion colli- 
sions. (3) LUCIAE introduces the suppression of strange 
quarks and effective string tensor to make some parame- 
ters related to product of strange quarks and string ten- 
sor in the JETSET depend on incident energy, size of 
system, centrality etc. Generally the LUCIAE can ex- 
plain many data very well [20,21]. 

First, we choose the particle multiplicity produced in 
each event as stochastic observable m and simulated 
10000 events for p+p at every energy point to obtain nor- 
malized multiplicity distributions correspondingly. Then 
we performed A-scaling. As a result, we find that 
the multiplicity distributions are approximately satisfied 
with the A = 1 scaling law, as shown in Figure 1. 

We calculated (m) = m iP( m i) an d ttirms = 
y/^2i m i {i f n*) 2 )P{mi) in terms of these normalized 
multiplicity distributions directly. While we fitted these 
distributions with Gaussian distributions to obtain \x and 
a. By comparing n with (m) and a with wirmSi we found 
that these distributions are basically Gaussian. Table 
1 shows the fitting parameters and information entropy 



{Hdirect or H gauss , which was calculated directly or by 
Gaussian parameter, respectively (see details in Ref. [22] 
and in the following). 

In order to better investigate the A-scaling, we defined 
a coefficient 
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which characterizes the validity of the A-scaling. We in- 
vestigat its dependence of incident energy with different 
A values. Figure 2 shows that the L-dependence on inci- 
dent energy from A = 0.5toA = 1.6. As a result, we find 




FIG. 1: A = 1/2 (a) and A = 1 (b) scalings for the multiplicity 
distribution of p + p reactions for the different beam energies. 
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TABLE I: Parameters in the direct calculation and fits for p 
+ P- 



that L is nearly a constant in whole investigated range 
(20-200 AGeV) when A=1.35, i.e., the system obeys the 
A-scaling law most suitably with A = 1.35 (see Fig. 3a). 

We also deal with another two systems of C+C (5000 
events) and Pb+Pb (1000 events). Similarly, their mul- 
tiplicity distributions basically obey Gaussian distribu- 
tions. In the same way, the best A-scaling is obtained 
with A = 1.00 for C+C and A = 0.80 for Pb + Pb, 
respectively. Figure 3b and 3c illustrate these results. 

Through A-scaling for different systems, we found that 
the distributions of particle multiplicity obey Gaussian 
distributions approximately, especially for C + C and Pb 
+ Pb systems and they obey the A-scaling law in a wide 
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FIG. 2: L as a function of energies. 



FIG. 4: Energy dependence of information entropy for p 
+ C and Pb + Pb, respectively. 



p, C 




p b + p b siKp^^^iiff 3 




04 >onM M 1 

1 1 1 1 III 1 





FIG. 3: A-scaling for p + p with A = 1.35 (upper panel), for C + 
C with A = 1.00 (middle pannel) and for Pb+Pb with A = 0.80 
(lower pannel) for the different beam energies. 



energy range for a given system, which indicates that no 
phase transition and no change of reaction mechanism 
exist. This is expected for the simulated data because the 
underlying particle production dynamics in the LUCIAE 
is a smooth function of beam energy without a phase 
transition. 

We calculate the respective information entropy in 
terms of these distributions of particle multiplicity with 
the method proposed by Ma in Ref. [22] . Figure 4 shows 
the dependences of the information entropy on incident 
energy for the p + p, C + C and Pb + Pb systems. The 
information entropy increases with the incident energy 
and with the sizes of the system monotonously. Also, 
the calculated values of Hdirect are consistent with the 
value H gauss obtained from Eq. (6) in the Gaussian dis- 
tribution. Again, no indication of phase transition exists. 

In addition, we use the distribution of strange particles 
multiplicity [23] and the distribution of the number of bi- 
nary collisions [24] as two stochastic observables to make 
the A-scaling plots and extract the information entropy 
in the same way. The similar A-scaling was obtained and 
a monotonous information entropy was also observed. 

In summary, we have demonstrated, for the first time 
to our knowledge, the A-scaling of charged particle mul- 
tiplicity, strange particle multiplicity and the number of 
binary collisions using simulated p + p, C + C and Pb 
+ Pb collisions from Ei ab = 20 to E lah = 200 AGeV. 
The LUCIAE 3.0, which includes the final state interac- 
tions, was used for the simulation. The A-scaling means 
that these observables obey a certain kind of universal 
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laws, regardless of beam energy and collision system. We 
found that the scaling values A for charged particle mul- 
tiplicity distributions are 1.35, 1.00 and 0.80 for p + p, 
C + C and Pb + Pb collisions, respectively. Moreover, 
the information entropy calculated from charged multi- 
plicity distributions increases with the beam energy and 
with the colliding system size monotonously. Both the 
A-scaling and the entropy values show no dis-continuity 
as a function of beam energy as expected because the 
LUCIAE has no change of particle production dynamics, 



while they are associated with a phase transition in the 
simulated data. We expect that the A-scaling and the 
entropy variable can be a valuable tool to search for pos- 
sible dis-continuities in nucleus- nucleus collisions associ- 
ated with the onset of a QCD phase transition. Further 
checks for the models which incorporate the QGP phase 
transition are in progress. 
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